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SUMMARY 

Several calculations that have aided operation of the Plum Brook Reactor are de- 
scribed. First, a simple mathematical model and the resulting calculations which were 

done to study fuel utilization are described. A decision to buy 240-gram uranium 235 
235 

(U ) fuel elements and to use them three times each was made on the basis of the study. 
*The power distribution in the mixed loading is discussed and is observed to be more uni- 
form than that in a uniform fuel loading. Charge life weight factors are then derived; the 
^derivation differs in approach and in results from previous derivations by others. Finally, 
a "xenon clock" on-line computer program in use at the Reactor is discussed, and an 
approximation useful for hand calculations of complex histories is derived. 


INTRODUCTION 

In regular operation of a high power test reactor like the Plum Brook Reactor, it is 
necessary to have efficient techniques of predicting the reactor charge life and power dis- 
tribution, the xenon transients after shutdown from high power, and the effects of pro- 
jected changes in fuel-element loadings or arrangements. Because of the short cycle time 
of the core (about 2 weeks) and the small amount of manpower available to do calculations 
for the reactor, the techniques must not be time consuming. It is helpful if the techniques 
can be set up for calculation on a small digital computer. 

Some of the methods devised at Plum Brook for predicting reactor performance are 
different in some respects from those used elsewhere. This report describes several of 
the methods: a model for an analysis of reactor fuel utilization, a derivation of charge 
life weight factors, and a "xenon clock" program for an on-line digital computer. Some 
of the results are presented and their significance is discussed. 



CORE DESCRIPTION 


The core consists of a 3 by 9 array of MTR-type fuel elements, cooled and moderated 
by water, with beryllium reflectors. The core volume (metal and water) is about 100 
liters and the core power is 60 megawatts (thermal). The core has 22 fixed elements with 
18 fuel plates each and 5 shim control rods with fueled followers having 14 fuel plates 
each. Control is with the shim rods in a bank. The normal method of operation is with 
new fuel elements placed at the core ends at startup and moved toward the center at the 
end of each reactor operating cycle. Thus operation is with a "mixed loading" of new 
and used elements at the start of each cycle. 


FUEL UTILIZATION STUDY 

A study was made to determine what gains in reactor charge life and fuel utilization 

could be made by using fuel elements of different loadings than the 200-gram uranium 235 
235 

(U ) elements being used at the time. The quantities of interest are the following: 

(1) Reactor charge life, which is the maximum amount of time, in mega watt -days (MWD)' 

that the reactor will operate with a given loading. Charge lives of the order of 900 mega- 
watt-days or greater were wanted. # 

(2) Critical height, which is the indicated control rod position at startup. This must 
remain above 15 inches (about halfway inserted) to avoid excessive flux peaking below 
the rods and to provide sufficient shutdown margin at startup. 

(3) Fuel utilization, which is the average amount of U fuel in grams used from each ele- 
ment before it is discarded. (The cost of an element is nearly independent of its new 

fuel loading, so that the average gram burnup per element is a good measure of the fuel 
utilization. ) The value had been approximately 80 grams per element with the existing 
mixed loadings of new and used 200-gram elements. 

(4) The maximum burnup which can safely be achieved. On this type of 
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element, it was about 144 grams of U (private communication with Mr. T. P. Hamrick 
of ORNL). This value, which was taken as a limit for this study, established an upper 
limit on the fuel utilization. Also, the maximum burnup per element differs the least 
from the average when elements are used the same number of times each. Therefore, 
the average number of times that elements are used should be an integer. 

(5) Flux perturbation on existing fluxes in experiment facilities was to be minimized 
and a flat power distribution maintained. This would be true if the relative fuel distri- 
bution at startup remained nearly the same. 

To study the variation of all these quantities as functions of fuel loading, the follow- 
ing equations were used: 
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With rods out, not including xenon or samarium, 


K 


eff 


- K ex = K *° P NL 


(1) 


With rods inserted (critical) including xenon and samarium, 


K eff 1 K ex AK rods + AK xe+sm 


Substituting equation ( 1 ) into ( 2 ) yields 


^rods + ^xe+sml'^NL - 1 


( 2 ) 


(3) 


where = rf. eP. For PBR, eP - 1. 0 and is not sensitive to fuel loading. Also, 

P NL is nearly constant for the fuel loadings of interest. From reference 1 K gx was 

measured to have a value of 1. 137 for a uniform loading of 168-gram elements. When 

26 

values of 77 = 2. 07, cr = 630 barns, and a total metal -water cross section of 1700 centi- 
,meters squared are assumed, equation ( 1 ) gives P NL = °* 6929 for that core loading. 

This is assumed to vary only with experiment reactivity worth for the 3 by 9 core. Thus, 




NL 


= 0.6929 (1 + P exp ) 


(4) 


For this study the value of P NL was 0. 6896. 

The shim rod worth varies with fuel loading. Most of the worth is due to the fast 
neutron thermalization inside the control rod water passages. According to reference 2, 
the total reactivity worth a' the control rods in the PBR with 168-gram elements (4365 g 
total) is 0. 362 AK/K, of wrich 0. 227 AK/K is due to fast neutrons and 0. 135 AK/K is 
due to thermal neutrons from outside the rods. The thermal flux is inversely proportion- 
al to fuel loading . Thus, :he relative rod worth R is 

0. 227 + 0. 135 r 3 -- 

r = V_M_ 

0.227 + 0. 135 


and 


AK r„ds = aK Jods R 


( 6 ) 
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r ref 


where AK~, qc j s was measured in the reference loading of 168-gram elements. 
In general, with rods out, not including xenon or samarium, 
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t/core 


/ 

J core 


S 25 (r)<p(r) dr 


S 25 (r)<p(r) dr + 


/ 

t/core 


S MW (r)«p(r) dr + 


/ 

J core 


(7) 


Z FP (r)<p(r) dr 



( 8 ) 


( 9 ) 


(10) 


4 



The value of F was determined to be 1. 192 by a criticality measurement on a depleted 

i FP 

core with rods out. This value of F is assumed to be constant, and a = 50 barns per 

RO 

fission. If these values are used, is computed and plotted in figure 1. 




t 



Figure 1. - Infinite multiplication with rods out against core loading for several element 
weights. No samarium or xenon. 

For rods out (end of cycle), 
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(12) 


1 + 0. 0247 K 


RO 60 


vs 


= 0. 6896 K 


RO 


core 


(13) 


K 


RO 

OO 


1. 505 + 


90 


VS 


= 1.52 


core 


(14) 


Operation will stop when K R ^ is reduced to 1. 52. This establishes the end of cycle 

OO 

value of core loading for a given fuel-element weight (see fig. 1). The charge life 
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is then T = (M - M^)/ 1. 27, where 1. 27 is the number of grams of depleted in 1MWD 
of operation. In figure 2, T is plotted against core loading for several element weights. 



Core loading at startup, M, g 

Figure 2. - Charge life as function of startup loading for several element weights. 

With rods inserted (at startup), 
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where 



S Sm 

V 

core 


and 



y 

N used Sm 

element 


(15) 


(16) 


(17) 


The Sgm per element equals about 16 centimeters squared for each used element at 
startup. The problem is to determine the number of used elements in the core at start - 
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up. It is noted that M - = 1. 27 T. There are 22 fixed elements with 18 fuel plates 

each and 5 control rod fueled followers with 14 fuel plates each. This makes 22+(14/18) 
5 = 25. 89 equivalent fixed elements in the core. Thus, the depletion per element in a 
cycle is 


1 27 T 

25.89 25. 89 


(18) 


If it is assumed that each cycle is the same as the previous cycle, then 


M = N^® w m 0 + N“ w (m 0 - Am) + N'“ w (m 0 - 2Am) + . . . 


new, 


(19) 


n-1 


M - M f 

M = nN”® w m 0 - N a ® w — - ^ I n = number of uses 


25.89 


I 

£ = o 


( 20 ) 


M = M n - (M - M-)/i + - + . . . + 

u 1 \n n n / 


( 21 ) 


; If each element is used three times, M = 1/2 (Mq + M^), and so on. Resulting values of 
N ei W are s ^ own * n li? 111,6 8. Since N^ e< ^ = 25. 89 - equations (15) to (17) can 

now be solved for AK^^ g and from figure 4 the starting critical height can be obtained 
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Figure 3. - Number of new elements at startup as function of core loading for several element 
weights. 7 




Indicated rod position, in. 

Figure 4. - Control rod calibration in clean core with 4365-gram loading. 

for each case. This condition is plotted in figure 5. Figures 2 and 3 are cross-plotted 
to obtain the average utilization in grams per element in figure 6. Then figures 1 to 6 
were cross-plotted to obtain estimates of core performance as a function of new element 
weight and number of elements added per cycle. An example is given in figure 7 for core 
performance with elements used three times each. 

Many observations can be made about the results. For example, figure 2 shows that 
the total core loading at end of cycle will always be about the same for given element 
weight and experiment reactivity worth. This has been observed experimentally. Also 
some of the variables behave in a nonobvious fashion, such as the average utilization in 
figure 6. The results indicated that operation with 240-gram elements used three times 
each would provide improved fuel -element utilization while also giving longer charge life 



Figure 5. - Critical height at startup as (unction of fuel loading 
(or several element weights. 
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without any significant penalities. This was verified by calculations using computer 
codes for two-dimensional diffusion theory. Therefore, the decision was made to use 
this type of operation. 


POWER DISTRIBUTIONS 

In addition to the greater fuel utilization, the mixed loadings have the advantage of a 
more uniform power distribution in the core. This effect is shown in figure 8, which 



Fuel-element position number 

Figure 8. - Comparison of power distributions in uniform and mixed loadings. • 

compares the power distribution among elements in the mixed loading of 200-gram ele- 
ments with that for a uniform loading. The maximum power density is reduced about * 
10 percent; this gives a greater heat-transfer margin so that the reactor can run at 
60 megawatts with control rods at deeper insertions than was the case with uniform 
loadings. 

A property of these loadings is that the distribution of thermal flux among elements 
is not sensitive to the loading or even to the control rod position. Evidently the increas- 
ed fission power at the ends of the core in the mixed loading is compensated for by the 
larger absorption cross section in the end elements. As a result, it is possible to cal- 
culate the power distribution (and the burnup distribution) among elements at the start of 
each cycle using 



F.M. 

l l 


E F.M. 
i l 


( 22 ) 


The use of this expression makes fuel accountability easy and aids in loading the reactor 
for a flat power distribution. 
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CHARGE LIFE WEIGHT FACTOR DERIVATION 


Before the start of each cycle, the charge life of the PBR is predicted using a method 
originally developed at MTR (ref. 3). The core to be used is compared to some previous 
"reference" core by the equation 


T = T ref + S w i<“i - M r ef,i> 


(23) 


The weight factors were derived differently than those of reference 3. The units of the 
weight factors were written as 


MWD 


MWD 


■ x 


Reactivity 


Gram Unit reactitity Gram 


(24) 


The perturbation expression for the latter quantity is 


Reactivity _ 
Gram 


/ 


Si'S. 


Gram 


f <Pf ( r )<P s ( r )| d V 
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a + 


Gram 


(pJr)cpJr) 


dV 


v f Js f ^(r)^ s (r) 
7core 


dV 


(25) 


A similar expression relating MWD to reactivity is 


Reactivity _ </core 
MWD 


r p^ f + 

J 

•/core 


<pf (r)<p s (r) dV - 


/ 

•/core 


1 6S 


a + 


dt 


<p s (r)(p s (r) 


dV 


v f |s f (pf (t)tp g (r) dV 
•^core L J 


(26) 


Dividing equation (25) by equation (26) gives 
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Gram 


W, = 


/ 


/ 6S_ 

# <p f + (r)<p g (r)dV a 


Gram 


1* 


(r)(pjr) dV 


(27) 


/ 

‘'core 


|6i/Z 


dt 


- <pj»<p s (r)| 


dV - 


/ 


(r)<p„(r) 


dV 


The change in the fission cross section due to a 1-gram change in the fuel density is 


6yS f „ N 0 °f 


Gram V e A 


(28) 


where V g is the volume of a fuel element. The change in absorption cross section per 
gram is 


6S. 


Gram U 


235 


i>5y' gyFP 6v 26 
1 + a v02j i + os a + os a 


v Gram Gram Gram 


(29) , 


This fission cross section at any time t is given by 


uZr (f) = i/S- (o) - 6i/S f (t) 

A i *•! 


(o) - 1 


6i/S* 


aM; 


25 AN? 5 

1 rp 1 


tGram U 


235 T 


AN 


25 
core. 


(30) 


f 25 


where AM? 1 * is the average number of grams depleted at i, AN? 5 /AN^ 5 re is the i 
fuel-element burnout relative to the core average burnout, and T is the charge life of 
cycle in megawatt-days. Differentiating equation (30) yields 


th 


di/S. 


f i 


-5 i/S* 


(AMf 5 


dt 


Gram U 


235 


AN ? 5 


AN 


25 
core; 


(31) 
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In equation (31), the quantity (AN? 3 /AN 23 ) equals P. and the quantity AM? 33 /T 

x cure x x 

equals the average number of grams depleted at i per total megawatt -days, which is 
(1. 27/27). The total cross section at i at any time after equilibrium xenon has been 
attained is 


.25 


MW 


Sa.(t) = saf a (t) + SaJ ViV¥ (t) + Saf e+Sm + Saf r + Sa 


FP 


.26 


(32) 


where 


Sa 25 (t) - Sa 25 (o) - i±J“ 6i/S f (t) 


sa: 


Xe+Sm 


(t) ' |^Xe +r l + y Pm] S f. (t) 


6S. 

i Oi 


f i 


u 


26 6S 'i 

Saf 6 (t) = — 

1 25 1 + a 


(33) 


Equation (32) is then differentiated to give dSa^(t)/ dt. Equations (28 to (33) were eval- 
uated with the parameter values shown as follows: 's: 


.235 


2.47 
24 


Neutrons per fission in U , v 

Avagadro’s number, Nq atom s/g -mole 0.6023x10' 

Fission cross section. Of CT a / (1 + a ) 

Atomic weight, A, g/g-mole 235 

Ratio of total capture to fission captures for U ,1 + a 1. 19 

Fuel-element volume, V , cm 0 10/27 

Rate of depletion of U 235 , g/MWD 1.27 


Yield of xenon. 


>Xe 


Yield of iodine, yj 

Yield of promethium, yp m 

Microscopic fission product cross section, a , b/fission 50 


0. 003 
0.061 
0. 0113 


.235 


694x0. 8862x0, 975 corrected 


Microscopic absorption cross section for U 

3 Maxv 

236 236 

Microscopic absorption cross section for U , a , b/atom 2.71 


a 235 b/atom to Maxwellian with plate disadvantage factor included 

' rkfir* on o 
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Then, they were substituted into equation (27) to give 


w i = 


2.07 


/-* 


<p_ dV - 0.935 


J fat 


dV 



(34) 


Two-dimensional multigroup diffusion theory calculations were run to obtain values 
of cpf, <pg, and <p s for use in equation (34). The core loading is shown in figure 9. The 
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Figure 9. - Fuel-element configuration and weight factor for 
mixed loading of new and used 200-gram fuel elements. 

calculations were done with all rods out and the core poisoned out to critical. The re- 
sulting weight factors are given in figure 9. These values differed significantly from 
those obtained by the method of reference 3, both in magnitude and distribution. The 
average of these values should be about 1. 27“^ megawatt-day per gram as shown in the 
FUEL UTILIZATION STUDY section; this is true for the values shown in figure 9. 
Therefore, these values were adopted for routine prediction of charge life of the PBR. 


POISON EQUATIONS 
"Xenon Clock" Computer Program 

The program calculates the poison transients occurring after shutdown in the PBR 
due to iodine-xenon and promethium-samarium buildup in the reactor core during before- 
shutdown operation. The program is divided into two parts. The first part samples the 
reactor flux and rod position every 10 minutes and updates the calculated iodine, xenon, 
promethium, and samarium concentrations while the reactor is operating. The second 
part, done on demand, predicts what the transient would be if the reactor were to scram. 
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In the first part, the equations are solved assuming that the flux is constant over 
the previous 10-minute interval. The new values of the concentration replace the pre- 
vious values and the computer reverts to other scan-log functions until the next compu- 
tation 10 minutes later . The equations used are 

-wv] 

1 - e I (p 


Th v y I S f 

1^2) - ® ^ 


Xe(tg) = e 


" 0f Xe^ t 2” t l^ 




“Xe " x l 


a 


Xe 


+ e 


-AjO^-tj.) 


I(t 1 )A I - Yj?f<P (y x + y X e) s f^ 


(35) 


a Xe ~ X I 


a 


Xe 


where 


C ’ "Xe = X Xe + °Xe^ 

* A similar set of equations describes the promethium -samarium chain. 

The second part of the program calculates the times after a scram at which the core 
would be just critical w’ith rods out. The reactor operator need just press a button on 
the console to have these times displayed for him. Figure 10 is a plot of the total poison 



Figure 10. - Calculated poison cross section as function of time after reactor scram for the Plum Brook Reactor. 
Operation, 3 days at 60 megawatts. 
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cross section as a function of time after shutdown from high power operation. The nota- 
tion S represents the maximum poison cross section which the rods can override if 
ex 

they are pulled clear out of the core. The problem is to calculate the times t ^ and tg 
shown in figure 10. Time t^ is the amount of time after scram which the reactor oper- 
ator has to restart the reactor before the poison transient builds beyond the control rod 
capabilities; tg is the waiting time before the reactor can be restarted if the operator 
fails to restart the reactor within t^. 

After scram the equations governing poison dynamics are 

-X T t 

I(t) = I 0 e 1 


Xe(t) = Xe Q e + 


X Xe t . X I I 0 L -X Xe t .~ X I l 


- e 


X I " x Xe 


(36) 


Sm(t) = SmQ + Ph1q( 1 - e 


" X Pm^ 


The xenon equation can be written 


Xe(t) = Xe^t) + I (jf 2 (t ) 


(37) 


and the samarium equation can be written 


Sm(t) = ShIq + Pmgf.j(t) 


(38) 


where 


f x (t) = e 


~ x Xe t 


- e > 

X I " x Xe 


(39) 


f,(t) = 1 - e 


” X Pm t 


Thus the after-scram poison transients are functions of the Xe, I, Pm and Sm concen- 
trations at scram and of the three time functions f fg, and f^. 

The total poison transient after scram is made up of the xenon and samarium tran- 
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sients. It is convenient to work in terms of the total poison macroscopic cross section 

S T (t) = £ Xe (t) + S Sm (t) = |a Xe Xe(t) + CT Sm Sm(t)J (40) 

The magnitude of the change in total poison macroscopic cross section after scram is 
ASrp = S T (t) - E t (o). However, the reactivity worth of the xenon is greater than if it 
were uniformly distributed in the core because the iodine concentration is greater in 
regions of greater -than-average worth. Measurements in the PBR have shown that the 
"effective” xenon- samarium cross section after a shutdown is more nearly 

! AS T (t) = K 1 S T (t) - K 2 S t (o) (41) 

where 

K x = 1.28 

i 

K 2 = 1. 02 

j % 

Expanding AErp(t) gives 

ASrp(t) = K^Sji(t) - K 2 2rji(0) = |c r X gXe(t) + crg^^Sm (tfj - K 2 (a Xg XeQ + ^gj^SniQ) 

I 

Let 


Xe(t) = Xe^t) + lQf 2 (t) 

(42) 

Sm(t) = SmQ + Pm^^t) 

(43) 


AS T (t) = K 1 (a Xe Xe Q )f 1 (t) + K^o^I^t) + K^a^Pm^t) + K^a^Sm^ 

" K 2^ a Xe Xe 0 + a Sm Sm 0> (44) 

The method of calculation is as follows: 

(1) The amount of poison the rods can override when drawn full out is compu- 

GX 

ted from figure 11. 

(2) The value of A2 T (t) is computed for time after scram increments from 0 to 
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.12 



Indicated rod position 


Figure 11. - Amount poison rods can override as function of rod position. 

12 hours. Time is calculated by locating the two values of A2rp(t) which bracket 
S ex an£ * by linearly interpolating between the time points which bracket t^. 

(3) Time tg is calculated as time t^ except that A2rp(t) is computed in incre- 
ments from t max to 100 hours. Values of f^, fg, and fg are included in table I. 


Approximation of Xenon Transient for Complex Histories 

An approximate hand- calculation method was devised for predicting the xenon poison 
transient after scram for complex power histories. The xenon equation after a scram is 


Xe(t) = Xe^t) + I^t) 
where fj and fg are given in equation (39) 

f ^o/^o ^o\ 

AXe(t) _ hV AXe 0 * W A1 0 _ aI 0 | * !q \ Xe 0 »0 / 

Xe(t) " f]l Xe 0 + f 2 I Q ~ I Q + Xe Q " 

t i~r ~ + 2 

A o 


(45) 


(46) 


for 60-megawatt operation (Iq > 10 Xe Q ), and 
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TABLE I. - TIME FUNCTIONS 


[Iodine decay probability, Xj, 0. 1037/hr; xenon decay probability, A Xe , 
0.07596/hr; promethium decay probability, Ap^, 0.01285/hr.] 


Time, 

t, 

hr 

-x„ t 
f j_(t) = e Xe 

x i L-W e -¥\ 

-AP t 
f 3 (t) = 1 - e m 

Mt) = \ 

V x Xe V > 

0 

1. 0000 

0.0 

0.0 

1 

.9268 

. 09479 

. 01276 

2 

. 8590 

. 1733 

. 02537 

3 

. 7962 

. 2376 

. 03781 

4 

. 7379 

. 2897 

. 05010 

5 

.6839 

. 3311 

. 06222 

6 

. 6339 

. 3633 

. 07420 

7 

. 5875 

. 3876 

. 08602 

8 

. 5446 

.4051 

. 09769 

9 

. 5047 

.4169 

. 1092 

10 

.4678 

.4236 

. 1205 

11 

.4336 

. 4263 

. 1318 

12 

.4019 

. 4254 

. 1428 

13 

. 3725 

. 4216 

. 1538 

14 

. 3452 

. 4153 

. 1646 

15 

. 3200 

.4072 

. 1753 

20 

.2188 

. 3484 

. 2266 

I 25 

. 1497 

. 2799 

. 2747 

30 

. 1024 

. 2162 

. 3198 

35 

. 07004 

. 1626 

. 3622 

40 

. 04791 

. 1200 

. 4019 

45 

. 03277 

. 08734 

. 4391 

50 

. 02241 

. 06286 

. 4740 

55 

. 01533 

. 04485 

. 5067 

60 

. 01048 

. 03178 

. 5374 

65 

. 007173 

. 02239 

. 5662 

70 

. 004906 

. 01571 

. 5932 

75 

. 003356 

. 01097 

.6185 

80 

. 002295 

. 007648 

.6422 

85 

. 001570 

. 005314 

.6645 

90 

. 001073 

. 003684 

.6854 

95 

. 0007345 

. 002549 

. 7049 

100 

. 0005024 

. 001761 

. 7233 
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Evaluating the integrals yields 


Kt s ) = 


r I S f 


topi + . . • + <p s b s ) 


(51) 


where the constants B n are dependent only on the time intervals. To find an equivalent 
constant flux which will produce the same I(t g ) over the time interval 0 — t 


Kt s ) - <pypf I 

Jo 


s e -Xj(t s -t') ^ = <py jSf ^ _ e -^t s j 


(52) 


Equating the two values of I(t g ) and then solving and assuming <p is proportional to core 
power P yields 


P = PjB'j + . . . - P n B^ + P^ 


(53) 


* where 


+4 


B 


B T =. 
n 


n 


_x #s ' V “ X I^s ■ Sl-l^ 
e - e 


1 - e 


" x I t s 


1 - e 


'^s 


(54) 


The coefficient B’ will be different for each power history. However, if the time inter - 
n 

vals are fixed beforehand and an average power over each interval is chosen, then the 
values of B^ are fixed and the equivalent average power over the 48 hours preceding the 
scram is 


P = PjB^ + . . . + P g B^ (55) 

The length of the intervals was chosen to give 

(1) Equal weight to each P n 

(2) Total 48 hours before scram 

(3) No fractions of hours 

The intervals and their weights are as follows: 
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Interval 

Weight, 

percent 

0-1 

9.92 

1-2 

8.91 

2-4 

15.31 

4-6 

12.49 

6—8 

10. 16 

8-12 

14.91 

12-16 

9. 87 

16-24 

10. 78 

24-48 

7.65 


Given these values and a set of calculations of AE Xe against t after scram from equi- 
librium xenon for neveral power levels, the xenon transient for any power history before 
scram can be calculated with an accuracy of 5 percent of the poison concentration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 29, 1967, 

120-27-04-54-22. 
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APPENDIX - SYMBOLS 


A 

atomic weight 

B n 

coefficients in xenon tran- 

sient approximation 


thermal flux at element i, 


N/ cm 2 - sec 

f 

thermal utilization, not in- 


cluding xenon and samar- 
ium 

H c 

critical height, in. 

I 

atom density of iodine 

i 

subscript referring to core 


position 

K eff 

effective multiplication 


factor 

K ex 

K eff with roc * s out 

k ro 

oo 

infinite multiplication with 


rods out, no xenon or 
samarium 

AK rods 

control rod worth 

AK ref 

AK rods 

rod worth in reference core 

AK Xe+Sm 

worth of xenon and samar- 


ium 

M 

core loading at startup, g 

M f 

core loading at end of 


cycle, g 

M i 

U 235 in element i, g 

Mf 

core loading of all new ele- 


ments, g 

m 0 

weight of a new element, 

„ tt 235 
g U 


Am average depletion of U 233 

in one fuel element in 
one cycle, g 

Nq Avogadro's number 

number of new elements 
in core at startup 

number of used elements 
in core at startup 

N 23 atom density of U 233 in 

fuel element i 

n number of times each 

element is to be used 

P probability of escape from 

capture while slowing 
down 

P effective core power over 

period 0 to 48 hr before 
scram 

Pm atom density of prometh- 

ium 

P^ fraction of core power 

produced by element i 

P NL nonleakage probability 

R relative rod worth 

Sm atom density of samarium 

T charge life of cycle, MWD 

t time, MWD or sec 

tmax time after scram at which 

poison concentration 
reaches a maximum, hr 

t time of scram 

s 
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u 

average utilization, g/ element 

o 

S MW 

macroscopic absorption cross 

V 

core volume, cm 0 


section of metal and water, 

V e 

O 

fuel-element volume, cm 


cm” 

w i 

charge life weight factor at 

2 25 

macroscopic absorption cross 
section of U , cm 


position i, MWD/g U 222 


Xe 

atom density of xenon 

r eq 

L Sm 

equilibrium cross section of 

i 




samarium, cm 

a 

ratio of captures to fissions in 

y eq 

^Xe 



y 235 

equilibrium cross section of 


yield , nuclei/fission 


xenon, cm’ 

7 


7 



A2 t 

change in total poison cross 

e 

fast fission factor 


_ 1 



section after scram, cm 

V 

neutrons/ capture in U 222 

a FP 

microscopic fission product 

X 

decay probability, sec”'*' 


cross section, cm 2 /g U 22 ^ 

V 

neutrons/fission in U 22 ^ 

°f 

microscopic fission cross sec 

k 

dummy variable 

tion, b 

p 

reactivity 

a 25 

microscopic absorption cross 
section for U 222 

p exp 

reactivity worth of experi- 



ments 

<P- 1 

above -thermal adjoint flux 

2 a 

absorption cross section, cm'* 

< 

thermal adjoint flux 

T 

core 

core cross section, cm”* 

<p(r) 

thermal flux at r, n/cm -sec 

S ex 

amount core can override 

^s 

thermal flux 

S f 

core fission cross section, 


2 

average thermal flux, n/cm - 


- 1 
cm’ 


sec 


2_ peak samarium after scram, 

s -1 
cm 

FP 

2 macroscopic absorption cross 

section of long lived low 
cross section fission pro- 
ducts, cm’^ 
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